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Abstract 

We use fermion mean field theory to study possible plaquette ordering in the 
antiferromagnetic SU(4) Heisenberg model. We find the ground state for both 
the square and triangular lattices to be the disconnected plaquette state. Our 
mean field theory gives a first order transition for plaquette ordering for the 
triangular lattice. Our results suggest a large number of low lying states. 
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Recently, the importance of orbital degeneracy in the physics of the insulating phase of 
transition metal oxides has been emphasized [0-0]. The understanding of the ground state 
and elementary excitations in these materials requires consideration of both the orbital and 
spin degrees of freedom. In V2O3 for example, the anomalous ordering of the antiferromag- 
netic insulator can be understood quite simply as a result of spin-orbital coupling 
The interplay between spins and orbitals may also be responsible for the lack of magnetic 
ordering down to very low temperature in L1N1O2 [PHl^|. The orbital ordering and corre- 



lations have recently been observed in synchrotron x-ray diffraction measurements []y| and 
in resonant x-ray scattering experiments [[14] . 

It has been pointed out that for a S = 1/2 system with a double orbital degeneracy there 
is an ideal limit where the spin and orbital degrees of freedom can be unified into a higher 



SU(4) symmetry |15],[|. While this symmetry is usually not realized in the transition metal 
oxides due to Hund's rule and anisotropy in hopping amplitudes, the SU(4) symmetric point 
may still shed light on the physics of more realistic systems, especially if the deviation from 
the ideal limit is not too large. Using fermion mean field theory(MFT), Li et.al. || showed 
that the SU(4) ground state on a square lattice possesses no magnetic long range order. 
Possible candidates for the ground state include the commensurate flux phase (1/4 flux per 
plaquette), and the plaquette solid state. The commensurate flux phase has broken time 
reversal invariance which is troublesome (although time reversal invariance may be restored 
when projected into the physical Hilbert space of one fermion per site). Simple variational 
calculations on the square and triangular lattice support the lack of long ranged order, 
but also suggest the ground state should be a plaquette solid or a resonating plaquette 
liquid H]. On the other hand, Schwinger boson mean field theory predicts a long range 



ordered state More recently, finite size numerical calculations on the square lattice 



provides further evidence that the ground state has no magnetic long range order [|i~7 
However, while the commensurate flux phase is two-fold degenerate, their numerical results 
suggest a 4-fold degeneracy. 

A possible explanation for the 4-fold degeneracy is the spontaneous formation of pla- 



quette state with alternating plaquettes of strong and weak correlations. This would be 
consistent with variational calculations. In this paper, we will investigate plaquette ordering 
within fermion MFT, a possibility not considered in the mean field calculation of Li et.al. [[|. 
Using fermion MFT, we will show that the uniform bond amplitude mean field ground state 
is unstable with respect to plaquette formation on square and triangular lattices. Plaquette 
ordering breaks the lattice translational invariance. Since the symmetry is discrete, a finite 
temperature phase transition is possible even in 2 dimensions. Within the MFT, we find 
the plaquette transition to be first order for the triangular lattice. For the square lattice, 
the result turns out to be ambiguous. 

We start with the antiferromagnetic SU(4) Hamiltonian j3|. 

m,n 

where (i, j) are the pairs of the nearest neighbor sites i and j on the lattice, and m = 1, 2, 3, 4 
is the flavor index. are SU(4) generators. The Hamiltonian is equivalent to an isotropic 
antiferromagnetic quantum spin-1/2 system with two-fold degenerate orbitals ||. In terms 
of electron operators, S^(i) = c im c ijn , where c ijH is the annihilation operator of an electron 
of flavor n at site i, the Hamiltonian in eq.(|l]) may be rewritten as 

H = Yl Yl 4m Ci ' n °\n C J,m- (2) 

m,n 

It is implied that there is a constraint on the number of electrons at every lattice site, given 

In this paper, we shall apply a fermion MFT to study the possibility of plaquette order- 
ing for the square and triangular lattices. Such an ordering would break the translational 



invariance of the crystal lattice. The fermion MFT is well known from literature flT8| , |I~9 
We define the bond operator §ij = J2 n c ln c j,m anc ^ introduce Lagrange multipliers to 
describe the constraint on site i. We then use mean field to approximate the bond operator 
by its self-consistently calculated average, and assume that (3>jj) = = Uj is re& l- Also, 
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we replace the local constraint by a global one and set Aj = A, independent of the site. The 
mean field Hamiltonian is then given by, 

,t 



H mf = -Yl ^ {in C h n + h.c) + 4 " A ( E C U C *> ~ N \^ ( 3 ) 

(i,j),n \ i,n 

where N is the total number of lattice sites, and an overall constant has been dropped for 
simplicity. We see that A acts as a chemical potential. 

We first discuss the square lattice case. As shown in Fig. ([]]), we consider two types of the 
bonds in the lattice, denoted by type- A and type-B, representing stronger and weaker bonds, 
respectively. In a square lattice, there are equal numbers of type-A and type-B bonds. The 
hopping amplitude on each bond is thus given by 

Uj = t a , if e {A} 

= t b , if (i,j)e{B} (4) 

The mean field Hamiltonian becomes 

H mf = - t ( cos k x a + cos k y a) cl n c %n 

k.ri 

+ irt (sin La) ct , cr 

Z-^l K X I fc+(7T,0),n «>« 

k,n 

+ ir t (sin k v a) ct , N ct„ 



k.n 

\ k,n ) 

In the above eqn., t = t a + t b , r = (t a — t b ) / (t a + t b ), k is a crystal momentum, and the sum 
is over the first Brillouin zone. 

In a square lattice, the Hamiltonian is invariant under each of the following transfor- 
mations: a) t — > -t, i.e. (t a ,t b ) -> (-t a ,-t b ); b) r -> -r, i.e. (t a ,t b ) — ► (t b ,t a ); and c) 
r — > 1/r, i.e. (t a , t b ) — > (i a , — t&). Therefore, it suffices to consider the parameter space t > 0, 
and < r < 1. The case of r = corresponds to the uniform bond amplitude mean- field 
state and the case of r = 1 corresponds to a periodic array of disconnected plaquettes. 



The mean field Hamiltonian can be diagonalized and we obtain 

Nt 2 I 



k,a,n \k,a,n 



where 01 is the creation operator for a particle in the state \k, a, n) with the energy e? 

k,a,n K,a,n 

Here, a = 1, 2, 3, 4 is the band index, and the sum over k is restricted to within the reduced 
Brillouin zone, which is a quarter of the first Brillouin zone of the crystal. The reduction of 
the Brillouin zone is the result of the translational symmetry broken by the plaquette state. 
The four flavor states for a given k and a are energetically degenerate. Numerically, we find 
that for r > 0.3, the four energy bands are separated by band gaps. The ground state of the 
system is a filled Fermi sea, with Fermi energy [if = X determined by the electron number 
equation, 



S^Wt t ) = V 6 ( A - e % ) 



N. 



k,a,n k,a,n 

To obtain the mean field ground state of the system, we minimize the energy E below 
with respect to t and r, 

- = -Y er efX-er )+-(l+r 2 ). (6) 

j\T j\T L. 4 k,a,n y k,a,n J 2 

k,a,n 

We first minimize E with respect to t for a given r. The energy as a function of r thus 
obtained is plotted in Fig.(Q). As shown in Fig.@, the uniform bond amplitude state, given 
by r = 0, i.e., t a = tj,, is locally unstable with respect to the plaquette formation. The energy 
decreases monotonically as r increases from to 1. The lowest energy state corresponds to 
r = 1, the disconnected plaquette state. The local stability of r = 1 state can be confirmed 
independently by solving the weakly connected plaquette problem using simple second order 
perturbation. 

The plaquette state is 4-fold degenerate. The discrete symmetry implies the ordering 
will remain stable at low temperature in 2D and vanishes either continuously or through a 
first order transition at some finite critical temperature. Within MFT, there will be another 



critical temperature above which t vanishes. Since t represents short-ranged order, this 
transition is unphysical and is an artifact of this type of MFT. 
At a finite temperature, the mean field free energy is given by 



The mean field amplitudes t and r are obtained by minimizing the free energy at a given 
temperature. As T is increased from zero, we find that the F vs. r curve is qualitatively 



1) remains as the lowest free energy state (note that since r = 1 is trivially a self-consistent 
solution, it is always a free energy extremum). Meanwhile the short-ranged correlation t 
decreases as T increases. In Fig.(H) we show the values of t that minimizes the free energy 
for r = and r = 1 as a function of T. Other r value curves lie between these two curves. 
We see that at critical temperature T t = 0.75, t —>■ for all values of r. Above T t , even 
nearest neighbor sites are uncorrelated. T t can be calculated analytically easily. The free 
energy Fmf for Hmf in eq.(|7]) can be formally expanded in a series in t^. Assuming the 
transition in t is continuous (which is likely since Fmf is even in t), T t is given by the 
temperature where the coefficient in front of the quadratic term changes sign from positive 
to negative. To quadratic order, Fmf can be easily found using perturbation by taking the 
chemical potential part of Hmf as the unperturbed Hamiltonian H and the remaining part 
as the perturbation if' (essentially doing high temperature expansion in H'). By expanding 
to second order in H', we obtain 



From this, we see that T t = 3/4 and that T t is independent of the configuration of t^, hence 
of r. The best configuration for ty is determined by the higher order terms in Fmf- Our 




where (3 = l/fc^T, and A is determined by the electron number eq., 




identical to the E vs. r curve at T = (Fig.(fJ)). Thus, the disconnected plaquette state (r 




(7) 
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calculation indicates it corresponds to r = 1. Since for all T < T t , r = 1, the transition in r 
is preempted by the transition in t. Within our MFT, we are therefore not able to answer 
the question of whether the plaquette ordering transition is a first order or a continuous 
one. However, the mean field transition in t is unphysical and merely signifies a crossover 
from weak to strong nearest neighbor correlations. We thus interpret our mean field result 
to imply that as soon as significant nearest neighbor correlations develop, the system has 
complete plaquette ordering. 

We now turn to the triangular lattice case. Similar to the square lattice, we consider 
plaquette ordering due to spontaneous formation of different (type-A and type-B) bond 
amplitudes as shown in Fig.(JJ). On the triangular lattice, there are more type-B bonds than 
A-type. Out of 3N bonds, 5N/4 are A-type, and 7N/4 are B-type. Rigorously speaking, 
the bonds within A-type (or B-type) are not all equivalent. For example, among the A-type 
bonds illustrated in Fig.(|l|), the diagonal one has different connectivity from the four edge 
bonds. This requires, in principle, to assign different amplitudes for these bonds in a fully 
self-consistent mean field theory. Instead, we approximate all A bonds and all B bonds to be 
identical among themselves, which simplifies the calculation greatly. Our calculation then 
amounts to a variational calculation with two types of bonds or equivalently, a mean field 
calculation that is self-consistent when averaged over the A bonds or over the B bonds. The 
physics should not be affected qualitatively by this simplification. Quantitatively, we have 
also checked this point by examining a single A-plaquette problem explicitly. Assuming all 
bonds to be equivalent, we find t a = 1.025, while allowing the 4 edge bonds (t' a ) and the 
diagonal bond ( t") to be different, we find t' a = 0.943 and i!' a = 1.333, giving an average 
of 1.021. The ground state energies in the two cases are also within 2% of each other, with 
values —1.31 and —1.33 respectively. The simplified scheme of only allowing two different 
bond strengths is therefore sufficient for our purpose. 

We thus introduce two average mean field amplitudes: one for A-type and one for B- 
type. The self-consistent condition is then given by their average one, t a = /Na, 
and tb = /Nb, where G {A} and G {£>} respectively, and Na and Nb are the 
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number of A and B type bonds. With this approximation, the mean field Hamiltonian in a 
triangular lattice reads, 

H m f = — 5^ ( cos ^c a + cos ky a ) + {1 — r/2)t cos (k x a + k y a)] c £ n c £ n 

k,n 

+ XT t ( Sil1 k * a + \ Sil1 + 4+(^,0),n C fc,™ 
k,n 

+ ir t (^sin k y a + i sin (k x a + k y a)^j 4 +(0i7r ), n c fc,n 

k,n 

— r t - cos (k x a + k v a) ct cr 

fc,n 

+ ^!(3 + 3^-r)-Afccl 3 ,„-iv]. (8) 

In the above eqn., we have used the same notation as for the square lattice. Note that in 
the triangular lattice, there is a lack of symmetry between the strong and the weak bonds. 
Consequently, the Hamiltonian does not possess all the symmetries of the square lattice. 
In particular, the system is not symmetric under the transformation r — > —r. From the 
definition of r and t, we have a symmetry in parameter space, r — > l/r,t — > rt. This allows 
us to study r within [-1, 1]. t —t symmetry is also not present for the triangular lattice. 
For r = and r — 1, the energy is minimized by t positive, and we assume this holds for 
other r 7^ also. The mean field Hamiltonian can be written in the diagonalized form, 

H mf =V t t pi (3 t + — [3(l+r 2 ) -r]-X f V M & — iV 1 . 

k,a,n \k,a,n / 

The ground state is given by a filled Fermi sea , whose total energy is given by 

I = ^ E W ( A " + J P (1 + r 2 ) ~ r}. (9) 

k,a,n 

Since there is no r — > — r symmetry, i?(r) 7^ E{—r). Our mean field result shows that 
£"(r) has lower energy for positive r. Similar to the square lattice, the ground state is the 
disconnected plaquette state corresponding to r = 1. In this case, the ground state is 
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12-fold degenerate. Again, the local stability of r = 1 can be confirmed by second order 
perturbation theory This also gives a heuristic understanding of why the energy is lower 
for positive vs. negative r. The r = — 1 state is a plaquette state such that neighboring 
plaquettes are connected by a single diagonal bond of equal strength as the bonds in the 
plaquette. Now imagine this connecting bond is not equal but much weaker, which will 
make this state quite similar to r less than but very close to 1. The local stability of r = 1 is 
thus consistent with it being lower in energy than r = — 1, and by continuity with positive 
r having the lower energy vs. negative r for at least a range of r close to 1. 
To study the finite temperature phase, we minimize the free energy 



As in the square lattice, r = 1 remains the stable state at low T (Fig.([|)). On the triangular 
lattice however, the transition in r is not preempted by the one in t, and we are able to 
study its phase transition. At T « 0.74, the curvature of the free energy curve at r = 
changes sign and r = changes from a local maximum to local minimum. When T is 
increased to T c = 0.782, r = and r = 1 states become degenerate in free energy and there 
is a first order transition from the disconnected plaquette state to the uniform state (or the 
other way round if we consider lowering the temperature). These results are summarized 
in Fig@. Since the lack of reflection symmetry in r allows for a cubic invariant in the 
Landau theory for its phase transition, the first order transition obtained is to be expected. 
Indeed, the behavior of F vs. r as a function of T is precisely what one would expect from 
such a Landau- Ginzburg free-energy functional. Upon further raising the temperature, the 
average nearest neighbor correlation continues to decrease, with the mean field transition 
to a completely uncorrelated state occurring at T t « 0.8. Unlike the square lattice case, 
on the triangular lattice, the transition in t occurs at different temperature for different r, 
i.e. T t = T t (r). This can be readily understood from the lack of t — > — t symmetry. As a 



N 



F 



N 



H 1 + e 



%, ,,r A )] + A + -[3(l+r 2 ) -r], 




(10) 
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result, the corresponding Fmf (eq.(|7|)) contains odd power terms, which in turn results in 
first order transition in tij. Naturally then, the transition temperature will occur above the 
second order temperature at 0.75 and will depend on the configuration of hence on r. As 
always, the transition in t should be interpreted physically as a crossover. 

In summary, within MFT, we find the uniform mean field state to be unstable with 
respect to plaquette ordering at T = for both the square and triangular lattices. For the 
triangular lattice, the plaquette ordering transition at finite temperature is first order, in 
agreement with Landau theory argument. For the square lattice, the plaquette to uniform 
state is preempted by a mean field transition in t, and MFT does not give an unambiguous 
answer to if the plaquette transition is first order or continuous. Taking the MFT on its face 
value, the transition is continuous as there is no entropy jump, although this is complicated 
by it being also a transition in t. Based on the reflection symmetry in r, Landau theory would 
also suggest a continuous transition. However, we cannot rule out a first order transition due 
to a negative quartic invariant in the Landau free energy functional or due to fluctuations. 
For the triangular lattice, the plaquette transition is first order as expected from symmetry. 

Within MFT, the ground states of the SU(4) model in both square and triangular lattice 
are found to be disconnected plaquette states (r = 1). Naively, this implies huge degeneracies 
in the ground states. Consider a lattice with periodic boundary condition. All the plaquette 
states with plaquettes covering all the lattice sites have the same lowest mean field energy. 
A simple lower bound for ground state degeneracy would then be D « 2^ for the square 
lattice, and D m S N ^ 12 for the triangular lattice. However, such plaquette states are not 
orthogonal to each other whether we consider the mean field states or the actual SU(4) 
states. They may not even be independent, so these lower bounds for degeneracy may 
not be valid, but should be taken as suggestive of huge degeneracy. Obviously, the exact 
SU(4) ground state will not exhibit saturated plaquette ordering (there will be for example 
resonance between plaquette states), and the actual ground state will only be 4- fold and 
12-fold degenerate respectively. Nevertheless, we may still interpret the mean field ground 
state to imply a very strong, although not saturated plaquette ordering, so that plaquettes 
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are only weakly coupled. We would then expect a huge number of low lying energy states 
in the energy spectrum. This is consistent with the preliminary results of Penc et al. [2t| on 
the square and triangular lattices using variational calculations on relatively large systems. 
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summer fellowship at University of Cincinnati. The authors wish to thank F. Mila and X. 
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FIG. 1. A schematic representation of plaquettes on a) square lattice, b) triangular lattice. 
Each solid line represents a strong bond connecting two lattice sites, with hopping amplitude t a . 
Each dashed line represents a weak bond with hopping amplitude tb- 
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FIG. 2. Energy per site, (E/N), in eqn.(||) is plotted as a function of r for a square lattice. 
E/N monotonically decreases as r increases from to 1. The lowest energy state corresponds to 
r = 1, the disconnected plaquette state. 
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FIG. 3. The optimal values of t as functions of temperature T for r = and r = 1 on a square 
lattice. For < r < 1, t v.s. T curves lie between these two curves. The short range correlation 
between neighboring lattice sites, characterized by t, decreases with increasing T. At T = 0.75, 
t — > for all values of r. 
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FIG. 4. Free energy per site from eqn.(|i~0|) is plotted with respect to r for the triangular lattice 
case at T = 0.1. At this temperature, r = 1 remains the stable state and r = remains as a local 



maximum. 
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FIG. 5. Free energy per site v.s. r on the triangular lattice, a) For 0.74 < T < 0.782, r = 1 
is the lowest free energy state while r = is locally stable. b) At T = 0.782, r = and r = 1 
become degenerate. At this T, the transition in t for a given r has already occurred for a range 
of r < 0. c) At T=0.783, r = is the stable state and r = 1 is a local minimum. Transition in 
t begins to occur for some r > 0. d) At T = 0.7836, the transition in t has occurred for r = 1, 
leaving r = as the only minimum. 
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